/*
 * MineralPhase.cpp
 *
 *  Created on: 24 May 2011
 *      Author: allan
 */

#include "MineralPhase.h"

MineralPhase::MineralPhase()
{}
	
MineralPhase::~MineralPhase()
{}
	
const vector<double> MineralPhase::ActivityCoefficients(double T, double P) const
{
	vector<double> activityCoeffs(GetNumOfSpecies(), 1.0);
	
	return activityCoeffs;
}

const vector<double> MineralPhase::Activities(double T, double P) const
{
	vector<double> activities(ActivityCoefficients(T, P));
	
	for(unsigned iSpecies = 0; iSpecies < GetNumOfSpecies(); ++iSpecies)
		activities[iSpecies] *= GetMolarFractionOfSpecies(iSpecies);
	
	return activities;
}
